Derivation of the Effective Chiral Lagrangian for Pseudoscalar Mesons from QCD 
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We formally derive the chiral Lagrangian for low lying pseudoscalar mesons from the first principles of 
QCD considering the contributions from the normal part of the theory without taking approximation. 
The derivation is based on the standard generating functional of QCD in the path integral formalism. 
The gluon-field integration is formally carried out by expressing the result in terms of physical Green's 
functions of the gluon. To integrate over the quark-field, we introduce a bilocal auxiliary field &(x, y) 
representing the mesons. We then develop a consistent way of extracting the local pseudoscalar degree 
of freedom U(x) in <&{x,y) and integrating out the rest degrees of freedom such that the complete 
pseudoscalar degree of freedom resides in U(x). With certain techniques, we work out the explicit 
[/(a:)-dependnce of the effective action up to the p 4 -terms in the momentum expansion, which leads 
to the desired chiral Lagrangian in which all the coefficients contributed from the normal part of the 
theory are expressed in terms of certain quark Green's functions in QCD. Together with the exsisting 
QCD formulae for the anomaly contributions, the present results leads to the complete effective chiral 
Lagrangian for pseudoscalar mesons. The final result can be regarded as the fundamental QCD definition 
of the coefficients in the chiral Lagrangian. The relation between the present QCD definition of the 
p 2 -order coefficient Fq and the well-known appoximate result given by Pagels and Stokar is discussed. 



PACS number(s): 12.39.Fe, 11.30.Rd, 12.38.Aw, 12.38.Lg, 



I. INTRODUCTION 

The study of low energy hadron physics in QCD is a 
long standing difficult problem due to its nonperturbative 
nature. For low lying pseudoscalar mesons, a widely used 
approach is the theory of chiral Lagrangian based on the 
momentum expansion and the consideration of the global 
symmetry of the system without dealing with the nonper- 
turbative dynamics of QCD W || . In the chiral Lagrangian 
approach, the coefficients in the Lagrangian are all un- 
known phenomenological parameters which are determined 
by experimental inputs. The number of the unknown pa- 
rameters increases rapidly with the increase of the precision 
in the momentum expansion. For example, the chiral La- 
grangian for pseudoscaler mesons with three flavors up to 
the p 4 -terms given by Gasser and Leutwyler || contains 
14 unknown coefficients. When the p 6 -terms are taken into 
account, there are 143 additional unknown coefficients ||. 



This kind of approach has also been applied to the elec- 
troweak theory [Q for studying the probe of the electroweak 
symmetry breaking mechanism [^[^). Since parity and CP 
are not conserved in the electroweak theory, there are even 
more unknown coefficients in the electroweak chiral La- 
grangian than in the case of QCD So far, this kind 
of study is of the level of finding out suitable processes 
at future high energy colliders to determine the unknown 
coefficients in the electroweak chiral Lagrangian and inves- 
tigating to what precision the determination can be. The 
relation between the coefficients in the electroweak chiral 
Lagrangian and the underlying model of the electroweak 
symmetry breaking mechanism is yet not known except for 
some very simple models JEJ . 

Further study on understanding the relation between the 
chiral Lagrangian coefficients and the underlying dynam- 
ical theory will be very helpful both in QCD and in the 
electroweak theory for reducing the number of independent 
unknown parameters which makes the theory more predic- 
tive. There are papers studying approximate formulae for 
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the chiral Lagrangian coefficients based on certain dynam- 
ical ansatz [Q, but the approach is not completely from 
the first principles of the underlying theory. Attempts to 
build closer relations between the chiral Lagrangian and 
the long distance piece of the underlying theory of QCD 
by considering the anomaly contributions with certain ap- 
proximation also exist However, several aspects of it 
imply that such kind of approach needs improvement, e.g. 
(a) the theory does not include spontaneous chiral sym- 
metry breaking, and the chiral symmetry breaking scale is 
put in by hand; (b) without putting in the chiral symme- 
try breaking scale, the obtained pion decay constant F n 
is proportional to an imposed very low (~ 320 MeV) mo- 
mentum cut-off on the underlying theory of QCD; (c) the 
positivity of F% depends on a careful choice of the regular- 
ization scheme. The approach in Ref. [[l0| does not contain 
the above problems. But in Ref. JTo| ], the approximation 
of large- N c limit is taken from the beginning and the ap- 
proximation of picking up only the local scalar and pseu- 
doscalar pieces of the color-singlet quark-antiquark bilocal 
operator arising from integrating the gluon-field is taken in 
the derivation. With the latter approximation, the formula 
for the p 2 -order coefficient Fq in Ref. [[[o| is expressed in 
terms of an imposed ultraviolet cut-off, and the formula can 
hardly be related to the well-known Pagels-Stokar formula 
for Fq lyj. Therefore, further improvement of studying the 
effective chiral Lagrangian from the fundamental principles 
of QCD is necessary. Actually, the study can be divided 
into two steps. The first step is to formally derive the ef- 
fective chiral Lagrangian from the fundamental principles 
of QCD and express the coefficients in terms of certain dy- 
namical quantities in QCD, which gives the QCD meanings 
of the coefficients. The second step is to calculate the re- 
lated dynamical quantities in QCD to obtain the values of 
the coefficients. This paper is mainly devoted to the first 
step. 

In this paper, we develop certain techniques with which 
we are able to formally derive the effective chiral La- 
grangian for pseudoscalar moesons from the first principles 
of QCD without taking approximation, and all the coeffi- 
cients are expressed in terms of certain Green's functions 
in QCD. Such expressions can be regarded as the funda- 
mental QCD definitions of the coefficients. As a simple 



example, we show that, under certain approximations, our 
p 2 -order coefficient Fq reduces to the well-known approxi- 



mate formula given by Pagels and Stokar II |. A systematic 
numerical calculation of the coefficients by solving the re- 
lated QCD Green's functions in certain approximation (the 
second step) will be presented in a separate paper [?]. 

This paper is organized as follows. Sec. II is on the fun- 
damental generating functional in QCD. We start from it 
and formally integrate out the heavy-quark and gluon fields 
to obtain a formal generating functional for the light quark 
fields. In Sec. Ill, we introduce a bilocal auxiliary field re- 
flecting the light meson degrees of freedom with which we 
can integrate out the light quark fields. Then we develop 
a technique for extracting the degree of freedom of the de- 
sired local field U{x) for the pseudoscalar mesons from the 
bilocal auxiliary field, and formally integrate out the the 
remaining degrees of freedom of the bilocal auxiliary field 
to obtain a generating functional for the local field U(x). 
In Sec. IV, we develop certain techniques to work out the 
complete U (x)-dependence of the effective Lagrangian in 
the sense of momentum expansion, and obtain the effec- 
tive chiral Lagranian which is of the form given by Gasser 
and Leutwyler In this process we obtain the QCD ex- 
pressions for all the coefficients in the effective chiral La- 
grangian. A discussion on the relation between the present 
QCD definition of the 0(p 2 ) coefficients Fq and the well- 
known Pagels-Stokar formula (an approximate result) |fl| 
will be given in Sec. V. Sec. VI is a concluding remark. 



II. THE GENERATING FUNCTIONAL 

Consider a QCD-type gauge theory with SU(N C ) local 
gauge symmetry. Let Ai (i = 1, 2, • • • , V, 2 — 1) be the gauge 
field, ipa V an d be , respectively, light and heavy fermion 
fields with color index a (a = 1, 2, • • • , N c ), Lorentz spinor 
index 77, light flavor index a (a = I, 2, • • ■ , Nf) and heavy 
flavor index a (a = I, 2, • • ■ , Vj-). For convenience, we sim- 
ply call ip® 11 the "light quark-field" , \f r £,*' the "heavy quark- 
field" and A 1 ^ the " gluon-field" . Let us introduce local ex- 
ternal sources J ap for the composite light quark operators 
^ <T, P , where a and p are short notations for the spinor and 
flavor indices. The external source J can be decomposed 
into scalar, pseudoscalar, vector and axial-vector parts 
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J(x) 



-s{x) + ip(x)7 5 + f (x) + $ (x)7 5 



(1) 



where s(x), p(x), v fJl (x) and a^(x) are hermitian matrices, 
and the light quark masses have been absorbed into the 
definition of s(x). The vector and axial- vector sources ^ (x) 
and jb (x) are taken to be traceless. 

Since the contributions from the anomaly term to the 
effective chiral Lagrangian has already been studied in Ref. 
]io|]9|] , our aim in this paper is to study the complete normal 
part contributions. So, in this paper, we simply ignore 
the standard CP-violating term related to the anomaly by 
taking the 9- vacuum parameter 8 = 0. 

Following Gasser and Leutwyler we start from con- 
structing the following generating functional 



Z[J] = J ThpVipTWTWVAn 

x exp i / d 4 x{£{ip, A^) + rpJip} 



VipVip exp < i / d xip(i$ + J)^> 



V^V^VA^ApiAf,) 



x exp < i Id x 



Cqcd(A) - -[F l (A 



2C 



-glfAi + (i0 -M- gA)$ 



(2) 



where Cqcd(A) = —^A % A l ^ u is the gluon kinetic en- 
ergy term, M is the heavy quark mass matrix, I? = 
•tp^-^i/j are colored currents composed of light quark fields, 
— ■^[F l (A^)] 2 is the gauge-fixing term and Af(A m ) is the 
Fadeev-Popov determinant. 

Let us first consider the integration over V^'D^'DA fl for 
a given configuration of i\> and tp, i.e. the current serves 
as an external source in the integration over V^V^VA^. 
The result of such an integration can be formally written 
as 



TWVVVA^f^) exp \i / d 4 x C QCD (A) 



^[F\A,)f - gl?Al + M gA)V 



exp 

n=2 



d A X\ ■ ■ ■ <i 4 x., 



xGi\:x (xx,---, x n )ig (xi) ■ ■ ■ it: (*»)> 



(3) 



where G)} is the full n-point Green's function of the A) - 
field containing internal heavy quark lines and with given 



sources Tj*. For simplicity, the gluon field integration in 
this paper is limited to the topologically trivial sector. In- 
clusion of topologically non-trivial sectors only changes the 
intermediate results but not the final result (p5|). 

By Fierz reordering, we can further diagonalize the color 
indices of the light-quark operators, and get 

Gj*\-:k>i, . . ■ ,x^i\{ Xl ){^) a , 0ll ^\{xi)] ■ ■ ■ 

xfe(*»)(^)a„,J B 7^£(s»)] 



d x 1 ■ ■ ■ d x n g ^pl---p2 i x li X\, - ■ ■ , x n , x n ) 

where G^.'.'.'^™ (xi, x[, ■ ■ ■ , x n , x' n ) is a generalized Green's 
function containing 2n space-time points. Then (|2|) can be 
written as 

Z[J]= I ■DipV^expil I d i xip(i(fi+ J)ip + 



n=2 



d x\ ■ ■ ■ d x n d x l ■ ■ ■ d x n : 



^ '■■pre (Xl, Xi, ■ ■ ■ , X n , X„) 

III. THE AUXILIARY FIELDS 



1. The bilocal Auxiliary Field 

For integrating out the light quark fields if) and ip, we in- 
troduce a bilocal auxiliary field $( ar >)( b (x, x') by inserting 
into (||) the following constant 

V<$> S (ncQ^M (x, x') - C"(x)^ c (x')\ . (6) 



We see from (6) that the bilocal auxiliary field 
$(a??)(f<) ( X; 2;') embodies the bilocal composite operator 
(x)ipa 1 (x') which reflects the meson fields. Inserting (^) 
into (0) we get 



Z[J] 



VipViJjV<S> 
x,5^iV c $( a ")^)(x,x') - C"(x)V^ C (x') 
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oo - 

+N C / d4 xi ■ ■ ■ d 4 x n d 4 x[ ■ ■ ■ d 4 x' n 



= ir [J,*,n c 



^pi—Pn ( Xl > X l } " ' ' ' Xn > X n) 



x^^(x ll x' 1 )---<f a ^(x n ,x' n ) 



(7) 



The (^-function in (0) can be further expressed in the Fourier 
representation 

6\N c ®(x,x') - ii(x)ip(x') 

2^j-|-gi J d 4 xd 4, x'U(x,x')- (^N c ^(x,x') — ip(x)ip(x')j 



With this we can integrate out the tp and i/ 1 fields and get 
Z[J] 

= J P$Pnexpi| -iN c Tiln[i<? + J - U] 

OO „ 

+N C / dA xi ■ ■ ■ d 4 x n d i x' 1 ■ ■ ■ d 4 x' n 

n=2 

(-i) n {N c g 2 Y 



^Pl—Pn i Xl 1 X ll ' ' ' 1 Xn ' X n) 



x$^(x 1 ,x , 1 )---$^(x n ,x' n )\, 



(8) 



where Tr is the functional trace with respect to the space- 
time, spinor and flavor indices. 
Let us define the classical field IL 



IL = J 



uS 



(9) 



where S is the argument on the exponential in (||). Let 
To [J, 3>,II C ] be the effective action for II C with given J and 
<£>. II C satisfies]] 

6T [J,$,II C 



dli7{x,x') 
Then T [J, $, n c ] is explicitly 



= 0. 



(10) 



*In the conventional approach, one usually introduces an exter- 
nal source J coupling to the field II. With this, the right-hand- 
side of @ euqals to - J fTJ. Eq.@ corresponds to taking 
J = 0. Similarly, when taking J = 0, the effective action 
ro[J,'I , ,II c ] euqals to the generating functional Wo [«/,$, J\ \j=q 
for the connected Green's fuctions. This leads to the left-hand- 
side of (Hit). 



J Vn exp iN c | - iTr ]n[i$ + J — II] 
+ / d 4 xi • • • d 4 x n d 4z x , 1 • • • d 4 x' n 

n «* 



n=2' 

A | Pi---p n V 1' 1' i^u A n/ 



With these symbols, we can formally carry out the integra- 
tion over the fl-field in (||), and express the result by 

Z[J} = Jv<$> exp{ir [J,$,n c ]}, (12) 



2. Localization 

Since we are aiming at deriving the low energy effective 
chiral Lagrangian in which the light mesons are approx- 
imately described by local fields, we need to consistently 
extract the local field degree of freedom from the bilocal 
auxiliary field $( ar >)( b O (x, x'). The extraction should be 
consistent in the sense that the complete degree of freedom 
of the mesons resides in the local fields without leaving any 
in the coefficients in the chiral Lagrangian. Otherwise, it 
will affect the validity of the momentum expansion In 
this paper, we propose the following way of extraction, and 
we shall see in Sec. IV that it is really consistent. 

The auxiliary field <!> introduced in (^) has such a prop- 
erty which allows us to define the fields a and O' related to 
the scalar and pseudoscalar sectors of $ as 

(fiVfi' + O'VO't) 11 ^) = {l)^ (K){av \x,x) 

(fiW - ft'Wt) ^) = ( l5 ) Cv ^ b ^ a ^(x,x). (13) 

Here the cr-field repesented by a hermitian matrix describes 
the modular degree of freedom, and the il'-field represented 
by an unitary matrix describes the phase degree of freedom, 
i.e. 



«t+(s) = <r(x), n'Hx)Q'(x) = 1. 



(14) 



As usual, we can define U'(x) = ft' 2 (x) which contains a 
[7(1) factor such that detU'(x) — e" 5 ^, where the deter- 
minant is for the flavor matrix. The unitarity property 
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of U'(x) implies that d(x) is a real field. We can fur- 
ther extract out the U(l) factor and define a field U(x) as 



U'(x) 



'U(x). It is easy to see that &etU(x) = 1. 



Then we can define a new field Q and decompose U into 

U(x) = Q 2 (x), (15) 

which is the conventional decomposition in the literature. 
This U(x), as the desired representation of SU(Nf) R x 
SU(Nf)L, will be the nonlinear realization of the peu- 
doscalar meson fields in the chiral Lagrangian. Note that 
the way of introducing the U (x)-ficld is not unique but is up 
to a chiral rotation which does not affect the final effective 
chiral Lagrangian since the chiral Lagrangian is chirally in- 
variant. The fields a and $ are intermediate fields which 
will not appear in the final effective chiral Lagrangian. It 
is straightforward to subtract a from the two equations in 
L3h and using (114) to get 



e 1 N f ftf(x)lxi[Pa$ T (x,x)]tf(x) 
= e~^n(x)tri[P L ^ T (x,x)}n(x), 



(16) 



where P R and Pl are, resprectively, the projection opera- 
tors onto the right-handed and left-handed states, the su- 
perscript T stands for the functional transposition (trans- 
position of all indices including the space-time coordinates) , 
and we have expressed the result in terms of £1. Eq.(|l6|) 
builds up the relation between &(x,x) and U(x) [or 0(a;)]. 
Taking the determinant of ([ll]) we can express 'd(x) in terms 
of $ T (x, a;) as 



det 



det 



tn[P R <f> T (x,x)} 



t n [P L <S> T {x, 



(17) 



where tr/ is the trace with respect to the spinor index. 



Eqs.(13)-([L6|) describe our idea of localization. To realize 
this idea in the functional integration formalism, we need 
a technique to integrate in this information to the gener- 
ating functional (fjj]). For this purpose, we start from the 
following functional identity for an operator O satisfying 
det O = det (cf. Appendix for the proof) 

VU S(U^U-1) 5(detU -1) 

xT[0) SiQO^Q-n^OQ^) 
— const, (18) 



in which J T>U S(UW — l)<5(det[7 — 1) is an effective invari- 
ant integration measure and the function J-[0] is defined 
as 

— !-r=det0 [ Va 6(0^0 -^a)5(a-^). (19) 
With the special choice of 



0(x) = e *"/ tn[P R $ T (x,x)], 



(20) 



which satisfies det O — detC^, eq.(18) serves as the func- 
tional expression reflecting the relation ( fL6| ) . Inserting ( |l8| ) 
and (|2^) into the functional (|lj) and taking the Fourier 
representation of the 5-function 



v 



^ e ~iN c J dx H-(r2e> t n-n t on + ) 



(21) 



we get 



Z[J] = J V<I>VUVZ S(U f U - l)6(detU - 1) 
x exp|ir [J,$,n c ] + +iN c J d A x 

E(x) (e~ i ^f'tf(x)tr l [P R § T (x, x)]Q){x) 



xtr, 



with 



-e l N J n(x)tn[P L $> T (x,x)]n(x) 



(22) 



n 



{det[tr ; P R $ T (a;,a;)]det[tr i P L $ T (a;,a;)]}2 

x yp ( 7(5^(tr / P i? $ T )(tr ; P L $ T )-at ( T N 
xS(a — (7^ ) 



(23) 



In (|22|), the information about the relation ( |16| ) has been 
integrated in. 

Next, we deal with the functional integration over the 
<&-field. For this purpose, we define an effective action 
f[fi, J,S,$ c ,n c ] as 
g if [n,j,H,$ c ,n c ] 

= J £>$ exp|iL [J,$,n c ] +iT I [<f>] +iN c J d 4 x 



xtrt 



E(x)(e 1 N f &(x)tii[P a <$> T (x,x)}rt(x) 



-e N f n(x)tn[P L <f> T (x,x)}n(a 



(24) 
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in which the classical field $ c (x, x') is defined as 



ds e ff[u, j,s c ,$ c ,n c ] 



4V 



(25) 



where S stands for the argument in the exponential in ( |24| ) . 
$ c satisfies 



af[n, j,H,$ c ,n c ] 



= o. 



(26) 



<9<^ p (x,a;') 

With these symbols, we can formally carry out the J V<& 
integration in (p2[) and obtain 



Z[J] = / VUVE 5(WU - l)5{detU - 1) 

xcxp{if[n, j,s,$ c ,n c ]}. (27) 

Here we have formally integrated out all the degrees of free- 
dom in $(x, x') besides the extracted local degree of freedom 
U{x). This localization is different from those in the liter- 
ature @. 

Similar to the above procedures, we can formally in- 
tegrate out the S-field by introducing an effective action 
S ef f [U, J, E c , $ c , II C ] as follows 

e iS etf [u,J* e ,* c ,n c ] = J 2?Sexp*f[n,J,S,$ c ,n c ], (28) 

where the classical field E c is defined as 

„ / VEEexpif [CI, J,S,$ C ,II C ] 



/ VE exp*f[fi, J,H,$ c ,n c ] 



and satisfies 



dS eff [U,J,E c ,$ c ,Il c 
d(E c ) ab (x) 



0. 



(29) 



(30) 



Then the ^.-integration in ( p7j ) can be formally carried out, 
and we obtain 



Z[J] = J VU 5{U ] U - l)<5(detC/ - 1) 
x exp{iS eff [U, J,S c ,$ c ,n c ]}. 



(31) 



We see from ( |3l| ) that 5 e / / [£/, J, S c , $ c , II C ] is just the action 
for U with a given 

Using d|), (H), @, @, © and @, one can further- 
show the following important relation 



^Note that H c , "I> c and Il c are all functionals of U and J through 

§, © and (!)■ 



dJ a f{x) 



u fix 



J VE^p(x,x) exp»f[», J,S,$ c ,n e ] 
/ VE expzf[f>, J,S,$ C ,II C ] 



NM"{x,x). 



(32) 



In (|32|) , the symbol <f> c denotes the functional average of <f> c 
over the field E weighted by the action T[Cl, J, E, $ c , II e ]. 
Eq. is crucial in the derivation of the effective chiral 
Lagrangian. 



IV. THE EFFECTIVE LAGRANGIAN 

To derive the effective chiral Lagrangian, we need to ob- 
tain the U(x)- and J-dependence of S e ff[U, J, E c , <fr c , n c ]. 
Note that S e ff[U, J, E c ,<fr c ,H c ] depends on U and J not 
only explicitly in (|2^) but also implicitly via 3 C , $ c and 
n c through @, @ and (g). The remaining task of 
the derivation of the effective chiral Lagrangian is to 
work out explicitly the complete U- and J-dependence of 
S e ff[U, J,E C , $ C ,II C ]. The procedure is described as what 
follows. 

First we consider a chiral rotation 

Ja(x) = [Cl(x)P R + tf(x)P L ] [J(x) + if)] 

x[Cl{x)P R + Cl\x)P L ], 
$£(x, !/) = P(x)P R + Cl(x)P L ] <& T (x, y) 

x[Cl\y)P R + Cl(y)P L ], 
U n (x,y) = [Cl(x)P R + tf(x)P L ] U(x,y) 

x[Cl(y)P R + Cl^y)P L ]. (33) 

The present theory is symmetric under this transformation. 
Since the S-field is introduced in (||l]) and the operator 
ClO^Cl — CvOCl' is invariant under the chiral rotation, there 
is no need to introduce Eq. Furthermore, since det Cl = 1, 



we can easily see from (17) that $n(x) = $(x). The explicit 
dependence of S e f / [U, J, E c , <I> C , II C ] on U (x) comes from the 
explicit Ct(x) (ft T (x))-dependence of f [Cl, J, <6 C , II C ] in (|J) 
[cf. (^ ) and (Pq)]. After the chiral rotation, this term 
becomes 



+iN c / drxtrf 



, "(JO 



S(x) e N / tr^^x,^ 



-e^?tr,[ft*E(x,aO 
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which no longer depends on U (x) explicitly. Therefore, af- 
ter the chiral rotation, there is no explicit J7(a;)-dependence 
of S e ff[U, J, S c , $ c , II C ], i.e. the complete [/(x)-dependence 
resides implicitly in the rotated variables with the subscript 
f2. For instance, the effective actions Tq[J, $ c , II c ], r/[$], 
r[il, J, S, $ c , II C ] and S e ff[U, J, H c , $ c , II C ] can be written 
as 



T [J, ITJ = r [Jn, $fi, n nc ] + anomaly terms, 

r z [$] = r z [$ n ], 

f[n, j, s, * C1 n c ] = f [i, Jo, s, $ 0c , n nc ] 

+ anomaly terms, 
5 e// [[/, J,S c ,$ c ,n c ] = S ef f[l, Ju,E c ,^n c ,Uuc] 

+ anomaly terms. 

From ( pi| ) we see that 

e if[l,Jn,S,*ne,n ne ] 

-/Ma-p{ff„[A,*,,H 1 J + ir I M 

+iN c J d 4 xtr ;/ {S(x)[-zsin^^ 
+ 75 cos^]$ n (x,a;)} 



(34) 
(35) 

(36) 

(37) 



N 



f 



(38) 



where tr;y denotes the trace with respect to the spinor 
and flavor indices. The anomaly terms in ([34]), ( [361) and 
( |37| ) are all the same arising from the non-invariance of 
Trln[z<^ + J — II] under the chiral rotation. Note that the 
functional integration measure does not change under the 
chiral rotation, i.e. T>$T>H = T>$>qT)Hq since the Jacco- 
bians from $ — > $q and II — » IIr> cancel each other. We 
see that the [/(^-dependence is simplified after the chiral 
rotation. 

The second approach is the use of eq.(^2|). As we have 
mentioned in Sec. II that we ignore the irrelevant anomaly 
terms in this study. Then after the chiral rotation, eq.(|32|) 
becomes 



dJ° p (x) 



u fix, anomaly ignored 



(39) 

We see from d39|) that once the Jn-dependence of $n c is ex- 
plicitly known, one can integrate ( |39"|) over Jq and get the 
U (x)-dependence of S e ff [1, Jo, H c , <&nc, IIq c ] up to an irrel- 
evant integration constant independent of U(x) and J(x). 



From that we can derive the effective chiral Lagrangian 
and the expressions for its coefficients. There can be two 
ways of figuring out the Jn-dependence of $o c - One is to 
write down the dynamical equations for the intermediate 
fields S c , $n c and Hn c and solve them (usually this can 
be done only under certain approximations) to get the Jn- 
dependence of these intermediate fields. The other one is to 
track back to the original QCD expression for the chirally 
rotated generating functional (||) through (|32| ) and (|2q ) by 
reverting the procedures in Sec. Ill and Sec. II, which 
can lead to the fundamental QCD definitions of the chiral 
Lagrangian coefficients without taking approximation. We 
take the latter approach in this paper. Because of the 8- 
function s(n c ^^(x,x') - (x)i><£° (x')\ in ©, 



we can express 4*q c (x,y) as 



(40) 



/ V^V^V^V^VA^VE^{x)^(y)e l ^^^^- A '^ 



where 



S[ip, i>, *, A, E] 

+H(— i sin — + 75 cos — )]?/>}. 



(41) 



N f Nf 

In and all later equations in this paper, the symbol if) is 
used as a short notation for the chirally rotated quark field 
tpn^- In (pl|), V iljf-ipip] and <&' are the quantities defined 
m (||) and (0) 

expressed in terms of quark fields, i.e. 



n 



det{—tTi c [ip R (x)ip L (x)}} 



x det{ — tr; c [ip L (x) ip R (x)] } 

1 V n 



*As an integration variable, with or without the subscript Q 
makes no difference. Once the classical equation of motion is 
concerned, distinguishing the rotated ipn from the unrotated ip 
will be necessary. 
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x5(a - cr) 



det 



,2iO'(i) 



det 



(42) 



(43) 



where tr; c is the trace with respect to the spinor and color 
indices. J43| ) implies that the range of ^'(x) is [0, 7r). 

Note that instantons contribute to both p2] ) and ( |43"| ) 
[ p^[ . The [/a(1) violating field-configurations only cause 
nonvanishing ??' but do not contribute to Ti[^--iptp]. 



With (^0|)-(^3j), one can integrate 
souces and obtain 



L N c 

over the rotated 



i5 e //[l,Jf!,H c ,$nc,nnc 



anomaly ignored 



= j Vi/jVipV^V^VA^VE exp<{ ilM— ipip 



1 - 



+iJ<Px{C{ilj,$,y^,A ll ) 

■d' <&' 
+H (_ is i n _ +75C os_)]^} 



(44) 



For realistic QCD (iV/ = 3), cos(tf'/N f ) does not van- 
ish. We can then shift the integration variable S — ► 
2 — ipn / cos(d' /Nf) to cancel the pn-dependence in the 
pseudoscalar part of (fli]). After carrying out the integra- 
tion over S, we obtain 



iS B f /[i,Jo,3 c ,*ne,nn e ] 



anomaly ignored 



= / VipVt/iVVVVVA,, 



t) 



0' 



S[ ip a { - isin— +75 cos—)?/' 



x exp<{ «r z [— ^V] 



N f < 
0' 



+V# n + i nj5 ~ sn - pn tan — }ip} 

JSf 



(45) 



In (fig), there is no po-dependence in the pseudoscalar chan- 
nel, and the po-dependence appears in the scalar channel 
as the combination sq + tan — . 



Eq.(45) shows that 5 e //[l, Jn, , S c , $n c , n^c] is fie QCD 



generating functional for the rotated sources so+Po tan 

cmd w«t/i a special parity odd degree of freedoms 
—ill) -0 b sin jf- +ip 75"0 fc cos j^- frozen. After making a fur- 
ther C/a(1) rotation of the tp and ?/> fields, the angle j^- can 



be rotated away and the frozen degree of freedom becomes 
just the pesudoscalar degree of freedom {[} a ^ip b as it should 
be since this degree of freedom is already included in the in- 
tegrating in of the [/-field. The automatic occurance of this 
frozen degree of freedom in the present approach implies 
that our way of extracting the [/-field degree of freedom is 
really consistent, i.e. nothing of the pseudoscalar degree of 
freedom is left outside U. After the {7a(1) rotation, Tj and 
the Jacobian due to the rotation will give rise to an extra 
factor in the integrand, which is the compensation factor 
for the extraction of the [/-field degree of freedom. From 
the point of view of the auxiliary field this corresponds 
to the contributions from integrating out the degrees of 
freedom other than U, say the a and rf mesons. 

Now we are ready to explicitly work out the effective chi- 
ral Lagrangian to the p 2 - and p 4 -order. As is pointed out 
in Ref. the vector and axial- vector sources should be 
regarded as 0(p) and the scalar and pseudoscalar sources 
should be regarded as 0{p 2 ) in the momentum expansion. 

1. The p 2 -Terms 

We first consider the p 2 -order terms. To this order, the 
anomaly can be safely ignored. Expanding ( ^5| ) up to the 
order of p 2 , we obtain 



^e//[l, Jn, S c , <I>n c , Hoc] 

p 2 -order 

d i xtr f [F ab {x)s ab {x) + F' ab {x)p ab {x)] 
+ { d i xd i zG abcd (x,z)a^ ab (x)a^ cd (z), 



(46) 



[0 (x)^ b (x)} 



where 
F ab (x) = 

F' ab {x) = -( [^(aOVAx)] tan 



Gf: d {x,z) 



ti'{x) 
N f 



[ip a {x) ltJ , 1 ^ b {x)][^ C {z) lljlb ip d {z)\ 



(47) 



and the symbol (O) for an operator O appeared in ([47]) is 
defined as 



O 



(48) 
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where 



xSi^p a (-ism— +75 cos— 

For F ab (x), translational invariance and flavor conserva- 
tion [l6| leads to the conclusion that it is simply a space- 
time independent constant proportional to 5 ab . So that it 
can be written as 



F ab (x) = F 2 B 5 



ab 



where 



F o B o 



For F (x), parity conservation ]17[] leads to 



F' ab (x) = 0. 



(49) 



(50) 



(51) 



For G ab „ (x, z), translational invariance leads to the con- 
clusion that it can only depend on x — z. We can further 
expand this dependence in terms of 6(x — z) and its deriva- 
tives. To p 2 -order, the derivative terms do not contribute, 
and the only term left is 5(x — z) J d i zG ab 1 ^ d (x, z). The co- 
efficient J d i zG ab 1 ^ d (x 1 z) is again independent of the space- 
time coordinates due to translational invariance. Then 
Lorentz and flavor symmetries imply that J d i xG ab J :d is pro- 
portional to g^ v 5 ad 8 bc . There cannot be terms of the struc- 
ture S ab S cd since this term is to be multiplied by a^ ah a^ cd 1 
and <Zq is traceless. Therefore the only relevant part of 
x, z) is 



Qabcd 



Gf v cd (x,z) = 5(x - z)g,J ad 5 bc F 2 
+irrclcvant terms, 



(52) 



where 



4(iV) 



d 4 x[G^,' abba (0 1 x 



-Gz, 



fi' ,aabb 



(0,x)} 



8(N 2 



f 



1) 



d 4 x 



& a (oh^ b (om b (x) %ll5 r(x)) 



-^([/(0)7 M 75^ a (0)][/(x) 7/i75 ^( a ;)] 



[i> (0) 7 W(0)] )([^ (x)7 M 75^ a (x)] 



+ 7 ^([/(0)7 M 75^ a (0)])([/(x i 7., :,,''(..■)■ 



(53) 



Note that there is no term like tr/[i;^] in (|46[). The reason 
is that there exists a hidden symmetry sn — > h^s^h, po — > 
h^pnh, et^ — > h^a^h, and vu — > h^v^h + hHd^h in which 
the vector source transforms inhomogeneously. So that the 
vector source can only appear together with the derivative 
id 11 to form a covariant derivative, and a hidden symme- 
try covariant quadratic form of the covariant derivative can 
only be an antisymmetric tensor [cf. ((5^)] which does not 
contribute when multiplied by a symmetric coefficient of 
the type of (||). 

With ©, (|l|) and (H) the effective action © is then 



p 2 -order 



^e//[l, Jq., S c , 3>fic, IIn c ] 

r 

= Fq / d 4 xtr/[aQ + B sn] 
= F 2 J d A xtr f j[V[/t][ V/i [/] 
+ is [[/(s-ip) + [/t( s + ^)] 



(54) 



where is the covariant derivative related to the exter- 
nal sources defined in Ref. The integrand is just the 
p 2 -order chiral Lagrangian given by Gasser and Leutwyler 
in Ref. Now the coefficients Fq and Bq are defined in 
([33]) and (|5(]) and are expressed in terms of certain Green's 
functions of the quark fields. These can be regarded as the 
fundamental QCD definitions of F 2 and Bq- 

2. The p 4 -Terms 

The p -order terms can be worked out along the same 
line. The relevant terms for the normal part contributions 
(ignoring anomaly contributions) are 



3 4 -order, normal 



d xtrf 



-JC 2 (d^a^ - d u a^ l ){d ll a n ^ - d u a n ^) 
+^3[ a h} 2 + ^ian a n a n,^an,v + K, b a 2 a tv f [a 2 n ] 

+K, s SQtXf[sn\ + ICgpft + ICwpntr f \pn] 
+K,ns n a 2 n + ICusntTflan} - JC 13 ,V£ v Vn 4l v 



(55) 
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where the covariant derivative and the antisymmetric 
tensor V£ v are defined as 

dW n = dW n -iv&a» n +ia» n v£ 

V,T =d^ n - iv£v» n + itftift , (56) 

and the fifteen coefficients /Ci,---/Ci5 arc determined by 
the following integrations of the Green's functions 



d 4 x x»'x v ' 



= [(3^1 - g vv +9^9^) 

+2K 2 g^g»' u ']5 ad 8 bc + irrclavcnt terms 



-- / d 4 xd 4 yd 4 z 



x ^ (Oh" 1 75V 61 (0)] (ar) 7 w 75 ^ (x)} 
x $** (y)-f*<y 6 1> ba (y)} (^)7^7 5 V b4 (z)] 

0263 £0364 £0461 \}_( n tJ-1^2 „M3/14 



tf*/"*)^ + .g" lAI V 2AI4 /C 4 ] 

^0264^0463^0361 jl^gMlM2^A'3M4 



+5 aib 



ga 3 b 2 ^0264^046! ^^gfJ-lH3 gH2H4 



_|_ja 3 64^0462 ^0261 jl^gMiM3^M2A'4 



+s aib 



fia 4 b 2 ^0263^0361 J_^A«iM4gM2M3 



_|_ja 4 6 3 ja 3 6 2 ^0261 J_ ^gHiH4gH2H3 

+g ^ g W)lC 3 + g^g^Ki] 
^oi6 2 ja 2 6! ^a 3 6 4 ^a 4 6 3 [gMiM2 gM3/*4 
+ ( 5 MiM3 5 M2M4 + gfi^^Ms)^] 

+(5 aib3 <5 a3bl (5 a2b4 (5 a4b2 [g M1M3 g M2M4 2A^5 

+ ( 5 M1M2 5 M3M4 + jWWjWW)^] 
^,50164^0461^0263^0362 [g^lM4^M2M3 2^ 5 



+ ( g MlM2 5 M3P4 + gflfS^M*)/^] J 

+irrelavent terms, 



-y rf 4 x^[^(o)^(o)][^ c (x)^(x)] 



d 4 x( [^(0)^(0)] tan ^ [^(z)^)] 



x tan 



tf'(:r) 



/ /c 



1 1 ^Ay^o^o)] 

x [^ 2 (x)7^75 V" 2 (*)] [^ a3 (y)7M75^ b3 (y)] 



—JCn(S aib2 S a2bs S asbl + ( 5 a i t ' 3 (5 a302 (5 a20 i) 
+/Ci 2 (5 aibl (5 a2b3 (5 a3b2 + irrclavcnt terms, 



ACl3_ 576(^-1)/^ 



^>9iiv9^v - 2gnn>g V v') 



c 



a:" V (0)7"V (0) ^ (sch^M 



9„ _ rpHlH2H3lJ-i 1 9_ T^Ml/ i 2M3M4 

^y^l^4.iJ^2^3 ± A "T" ^!J^1^2 i/A I 3A l 4 J - A 

_„ „ rpHHJ,2H3l^4 _ 9 TiA'lM2M3A'4 

yp.HJ.3y^2P.i ± A ^yMl ^4 i/M2 M3 B 

9 ~ „ T^MlM2A'3A'4 1 TiA I lM2M3A'4 

^.y/Ji/J2. ( /M3M4- t _B ' y^i^yp.i^i-'- b 



K 



15 



4(7V2 - 1) 
x ; 6 (x)7 AI 75?/' a (x) 



V> (0)V> b (0) tan 
1 / 



t?'(0) 
7V f 



( V (0)^ a (0)tan 



i?'(0) 
(57) 



with Ta, Tg defined as 

d 4 xd 4 2 /x^^ ai (0) 7 ^?A bl (0) 

x^ 2 (x) 7 ^ 2 75^ b2 (x) V^ 3 (y)7 Al3 75^" 3 (») 
__ ^0162^0263^036121^1^2^3/^4 _|_ ^0163^0362^026121^1^2^3/^4 

+irrelavent terms. (58) 
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In ( B7p and (p8[), (■ ■ -)c denotes the connected part of {■ ■ •), 
and the irrelavent terms are those leading to tr/[a^] or 
try[u^] after multiplied by the corresponding sources. 

To further evaluate the rotated source parts in (|55|), we 
make use of the p 2 -order equation of motion 

1 



-B [p n - — tr/(p n )] 



(59) 



and the following identities 



■[n(s-ip)n + rt(s + ip)rt] 



Pa = - *p)n - n^* + ip)n f ] 

+ [-(V £/)[/*( Vt/) 



(60) 



where F^" and F£ are, respectively, the field-strength ten- 
sors of the right-handed and left-handed souces defined in 
Ref. §. With (||) and (fcj) and taking N f = 3, eq. (|f) 
becomes 



p 4 — order, Hormal 



J d 4 x L[ 


norm) [tr f (\7^\7 n U)} 2 


l (norm 




^ (norm 
+^3 


^[(V^V^) 2 ] 


(norm 




l (norm 


tr / [V^tV At (7(x t C/ + t/ t x)] 


. (norm 




(norm 


[tr/^CZ-X^)] 2 


^ (norm 


tr^xW^ + X^W] 


- t~ (nom 


'^/[F^VWt/t + F^ VELVET] 


j (norm 
+^10 








TT (norm 

+H 2 





where x = 2Fo(s + *p) 



(61) 

The integrand in (|6l]) is just 



the normal part contributions to the p 4 -order terms in the 
chiral Lagrangian in Ref. ||, and the coefficients are now 
defined by 



(norm 
2 

(norm 
3 

(norm 
4 

(norm 
5 

(norm 
6 

(norm 
7 

(norm 
8 

(norm 



_L/C 4 + i/C 5 + 1^3-^4, 
— (/C 4 + /C 6 ) + -/Ci3 - 7^^i4, 



2/C4 — 6/Cij 



16 

3/C14), 



16Bb' 
16B6' 

16JV) 



10 



£15 



16F 2 



16 



[ACi 



-(4/C 



13 



£14), 



16F VV/ 

r /Cg + 



1 „ I 



B, 



-(/C 2 - /C13), 
-7(^2 + £13), 



1 



(62) 



The twelve standard coefficients L\ 



(norm) 



L 



(norm) 



j (norm) jj(norm) 
L 10 ' U l ' 



jj(norm) expressed in terms of 



twelve independent p -order coefficients, IC2, /C3.4 = /C3 — 
2/C4, /C 4) 5 = /C 4 + 2/C5, /C 4) 6 = /C 4 + /C6, /C7, /Cg, ^1,9,15 = 
^1 — + jJ^15) /C14045 = X4 + "gr^lO + B^^Cl5; A^n, 

^Ci2, A^i3, and /Ci 4 . 

The total coefficients are then 

u = L (norm) + Wanom) f = X 1Q 



(63) 



where ]}-°- nom ^ anc j j^(°- nom ) are an0 maly contributions 
to the coefficients given in Ref. jl(],[)) . 

So we have formally derived the p 4 -order terms of the 
chiral Lagrangian from the fundamental principles of QCD 
without taking approximations and have expressed all the 
coefficients in terms of the integrations of certain Green's 
functions in QCD. Eqs. fiSTQ, $5$) and ^6^ ) give the fun- 
damental QCD definitions of the the twelve coefficients 



§In Ref. J2), x is defined as \ = 2Bq(s + ip)ez 9 . In this paper we have taken 8 — 0. 
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(norm) 



■L 



(norm) j^-(norrn) 



10 



and H. 



(norm) 



The procedure 



can be carried on order by order in the momentum expan- 
sion. 

The expressions @, (§|), @, @ and @ are conve- 
nient for lattice QCD calculations of the fifteen coefficients. 

V. ON THE COEFFICIENTS Fg AND B 

So far, we have given the formal QCD definitions of the 
fourteen coefficients of the chiral Lagrangian up to the p - 
order. To get the values of the coefficients, we need to solve 
the relevant Green's functions which is a hard task, and 
we shall present the calculations in a separate paper fl2| |. 
To have an idea of how our present formulae are related 
to other known approximate results, we take the p 2 -order 
coefficients Fq and B [eqs. (pOj) and (|5|)] as examples and 
make the following simple discussion. 

As we have mentioned in Sec. IV, there can be two ways 
of figuring out the explicit Jo-dependence of <&n c for evalu- 
ating S e f f from . One of them is to solve the dynamical 
equations for the intermediate fields $fi ci IIq c and S c , and 
the other is to track back to the original QCD generating 
functional without the intermediate fields. For convenience, 
we took the latter way in the above derivation of the chi- 
ral Lagrangian. To compare our results with the known 
approximate results, we are going to take certain approxi- 
mations, say the large N c limit, with which the calculation 
of the intermediate fields becomes even more convenient. 
Thus we take the former way in the following discussion. 

First we take the large N c limit. It can be easily checked 
that, in this limit, the functional integrations in (pd|), ( p^ ) 
and ( p8| ) can be simply carried out by the saddle point 
approximation (taking the classical orbit in the semiclas- 
sical approximation). The saddle point equations (|To|), 
( p6| ) and ( ^0| ) are just the dynamical equations determining 
Ilac'I'ncS,, as functions of Jo, which are 

^c ){b °^y) = +Jn--nn c )- 1 }W^(y,x), 



OO 



,4 .4 ,4 , ,4 , H)" +1 (iVcg 2 )" 

, a x\ ■ ■ ■ a x n d x l ■ ■ ■ a x n 



n=\ ' 

X(jr ppi-p„ \ x i Vi X li *1) "' i x ni x n)^Slc \ X 1> X l) ' 

x^ c p "(x n ,x' n ) = 0, 



tr;[(— isin 



&(x) 



- 75 cos 



tf(x) 



Nf N f 
where S is a short notation for the following quantity 

o r . My) 



(64) 



E">(x) = 



d^ c p (x,x) 



d j/ tr^/{S c (y)[-^sin- 



^/ 



+ 75 cos^]^ c (y, 2/ )} 



h c fixed 



(65) 



In @ and @, § c (x) depends on $(x,x) through (|l7|). 
Note that the effective action Tj [$n] belongs to 0(1/N C ), so 
that it does not contribute in the present approximation. In 
( p4[ ) , the field IIn c can be easily eliminated and the resulting 
equation is 

[i(j> + «$oc _1 + i> n + i f275 - s n + ipnj5 

+ ~} ap {x, y) + J2 ^ ' ' ' dix nd 4 x[ ■ ■ ■ d 4 x' n 

n=l J 

{-i) n+1 {N c g 2 ) n ^^-„ nt 
x^(x 1 ,x' 1 )---^(x n ,x' n ) = 0. 



X n) 



(66) 



In order to compare our results with the usual dynami- 
cal equations in the ladder approximation, we further take 
the ladder approximation which, in the present case, corre- 
sponds to ignoring all the n > 1 terms in ( |66[ ) and with 

G PiP2 {xi,x 1 ,x 2 ,x 2 ) 

= -\Gn 1 u 2 (xi,X 2 )(-f' 11 )a 1 p 2 (^ 2 )^ P l 

x5(x' 1 — X2)S(x' 2 — xi) + 0( — ) term, 

-Ac 

where G^ v {x,y) is the gluon propagator without internal 
light-quark lines. Then, in the ladder approximation, ( |66| ) 
becomes 

[i(j> + i^' 1 + i> q + 4 075 - sq + ipnj 5 
1 



+E](x,y) + ^NcG^ix^^x^Y = 



(67) 



On the other hand, in the large N c limit, $q c is just <f>oc 
which is the full physical propagator of the quark with the 
rotated sources. When the sources are turned off, <&n c can 
be expressed in terms of the quark self-energy E(— p 2 ) and 
the wave function renormalization Z(—p 2 ) by the standard 
expression 
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T(a»j)(60 







(68) 



[*SJ ( °" )( *W) 



«n=pn=e£=a&=0 



= 6 a6 



d 4 p 
(2tt) 4 



-ip(x-y) 



in which translational invariance and the flavor and parity 
conservations have been considered. Plugging (|S8|) into (|67[) 
we have 



xG llv (p-q)'f- 



1 



-7 V = 0, (69) 



-£(-g 2 ) 

where G liV {p) is the gluon propagator in the momentum 
representation, and the fact 5| sotirces= o = is taken into 
account. ( )69| ) is just the usual Schwinger-Dyson equation 
in the laddar approximation. 

With the solution of the Schwinger-Dyson equation, the 
formula ( |5Q| ) for FqB can be expressed as 

^■^/^ w^w (70) 

By definition [cf. (fig ) and (p2])1, Fq is related to the coeffi- 
cient of the term linear in an in the expansion of <&n c . We 
denote 



(71) 



[*L] (a " )(bC) (-,y) 



linear in ajj 



d 4 z 



d 4 pd 4 g 



Then Fq is determined by0 



(72) 



To this order, we still have S|j mear m qM = 0, and eq.(J67|) 
reduces to 



1 



d 4 fc 



/ T^G^Ak-qh^l^kW =0. 



(73) 



'"This is an alternative expression for Fq equivalent to (f 



In the literature, a further approximation of dropping 
the last term in ffi§J is usually taken |Q (It can be shown 
that to leading order in dynamical perturbation this 
term can be reasonably ignored |p^|.). Moreover, to leading 
order in dynamical perturbation or in the Landau gauge, 



Z(-P 2 



1. Then PJ? becomes 



F 2 



iNc 
8 

-UN, 



d A q 
d 4 q 



tr;[7 M 7 5 - 



1 d 



■7m75- 



(2tt) 4 [8dqK V? 2 -E 2 (-g 2 ) 



[£(-g 2 ) + ±g 2 £'(-g 2 )]£(-g 2 ) 



S(-g 2 ) J 



(74) 



[g 2 -£ 2 (-g 2 )] 2 

The anomaly contribution to Fq calculated in Ref. jl(| is of 
the same form as the first term in ( [74] ) but with an oposite 
sign, so that it just cancel the first term in in (f74|). Then 
( [74[ ) is just the well-known Pagels-Stokar formula for Fq 
p"T| . Thus the Pagels-Stokar formula is an approximate 
result of our formula by taking the approximations of the 
large N c limit, the ladder approximation and dropping the 
last term in j73] ) (or to leading order in dynamical pertur- 
bation). 



VI. CONCLUSIONS 

In this paper, we have derived the normal part contribu- 
tions to the chiral Lagrangian for pseudoscalar mesons up 
to the p 4 -terms from the fundamental principles of QCD 
without taking approximations. Together with the anomaly 
part contributions given in Ref. Jl^,^], it leads to the com- 
plete QCD theory of the chiral Lagrangian. 

We started, in Sec. II, from the fundamental generating 
functional (0) in QCD, and formally expressed the integra- 
tion over the gluon-ficld in terms of physical gluon Green's 
functions. Then we integrated out the quark-fields by intro- 
ducing a bilocal auxiliary field &(x,y) [cf. @]. To extract 
the degree of freedom of the local pseudoscalar-meson-field 
U(x), we developed, in Sec. Ill, a technique for extracting 
it from the bilocal auxiliary field $(2;, y) [cf. (|l3| ) and (jlH)] , 
and integrating in the extraction constraint into the gener- 
ating functional [cf. (^8|), (po|) and (^3)]. This procedure is 
consistent in the sense that the complete pseudoscalar me- 
son degree of freedom is converted into the [/(ar)-field such 
that the pseudoscalar degree of freedom in the quark sector 
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is automatically frozen in the path-integral formulation of 
the effective action S e ff. 

We then developed two techniques for working out the 
explicit [/(^-dependence of S e ff in Sec. IV. The first one 
is to introduce a chiral rotation (|3^) which simplifies the 
[/(^-dependence in such a way that the [/(cc)-dependence 
resides only implicitly in the rotated sources and some ro- 
tated intermediate-fields, and the second one is to imple- 
ment eq.(|3^) to obtain S e ff from the averaged field &n c . 
To avoid dealing with the intermediate-fields, we tracked 
back to the original QCD generating functional with which 
the implicit /7(x)-dependence only resides in the rotated 
sources. With all these, we expanded S e f / in power series of 
the rotated sources and explicitly derived the p 2 -terms and 
p 4 -terms of the chiral Lagrangian for pseudoscalar mesons 
[pj. In this formulation, all the fifteen coefficients in the 
chiral Lagrangian are expressed in terms of certain Green's 
functions in QCD [cf. @, (||), @ and @]. These for- 
mulae can be regarded as the fundamental QCD definitions 
of the fifteen coefficients in the chiral Lagrangian. These 
expressions are convenient for lattice QCD calculation of 
the fifteen coefficients. 

To see the relation between our QCD definition and the 
well-known approximate results in the literature, we took 
the p 2 -order coefficients Fq and Bo as examples in Sec. V. 
With the approximations of large N c limit, ladder approx- 
imation, and dropping the momentum-integration term in 
the approximate Bethe-Salpeter equation ( [73^ , our formula 
reduces to the well-known Pagels-Stokar formula for the 
pion decay constant fl]|j . The complete calculation of 
the fifteen coefficients will be presented in a separate pa- 



per |12|. The derivation of the effective chiral Lagrangian 
further including the n' meson or the p meson, and the ap- 
plication of the present approach to the electroweak chiral 
Lagrangian are all in preparation. 
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APPENDIX: 



Here we give the proof of eqs.(|18|) and (|19|) in the text. 
Consider a matrix operator O satisfying det O = det , 
we calculate the following functional integration 



I = J VU S(WU - l)5(detZ7 - 1) 

where T>U 8{WU — l)S(detU — 1) serves as the invariant 
integration measure at the present case. The two delta- 
functions 6(U'U—1), 5(detU— 1) constrain the integration 
to the subspace with unitary and unity determinant of the 
[/-field. 

We can rewrite S(detU — \)8{WU — 1) as 

5{detU - 1)5(UW - 1) 

= i<5([det[/] 2 - l)6>(detC/)(5([/ t [/ - 1) 



~<5(det?7[det?7 - det[/ t ])6»(det[/)(5([/ t [/ - 1) 

6(U*U- 1), 



Ldet[/-det[/t)^ detC/ ) 



detU 



so that 



1 = 2 I VU ~ 1 ^ detU ~ deiUt ) 

Next, we introduce two auxiliary fields E and E and write 
I as 

I = - J £>[/X>E£>E 5{U^U - \)8{detU - dett/ 1 ) 
x ^^p- S{t - E) F[0] <5(E - fitOfit) 

x<S(det[/ - detu rf^^ F[Q] 
x£(Q f [E - E]fi)<5(n f Eft - U^O) 

x8{9)tn - cftu). 



We then change the integration variables E and E into 

E -> E' = fttEft E -> E' = n f Efi, (Al) 
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and get 

I = \J VUVTlVf! 5(U^U - l)<5(detE7 - dettft) 

x 9(detU) ~, _ _ uf 

detU 
xS(t' -0+fJ) 

= -J VUVYlVYl 5(U^U - l)<5(detE7 - detrf) 

x||^(£'-£') T[0\ 5{UH'-0) 

x£(£'-E't) 
= ~ /" VUVYlVYl det(E' 5 )5(£'[[/+[/- 1]E') 

x<5(det[E'£/] - det[E'C/t])_g^)_j ( s' - £') 

x J"[0] <S([/E' - 0)5{t' - £*) 
= i /" P[£/£']PE'PE' det(E ,4 ) ( 5(E' t [/ t C/E' - E'+E') 

x<5(E' - E') JF[0] £([/£' - 0)£(E' - E't). 
We further change the integration variables {/ into 



U^U' = ITS', (A2) 



and get 

J =\J VU'VYlVY! S(U n U' - E'+E') 

x<5(E' - E') .F[0] 5(U' - 0)8(t' - E't). 
The [/' and E' integrations can be carried out and we obtain 

I = 2{detO}2 d{det ° ~ det ° f) / VY! 6{&<D ~ S ' tS,) 
xdet(£' 4 )0(g|^[0]«5(E'-E't) 

5(1 - / £>£' 5(0+0 _ s'tE') 



2[det0] 



detO 



:[det(E' t E')] 2 t 



detO 



=)Fp]8(Z' - E ;t ) 



v/det[E'tE']' 

_ ^!) / x>£' 5(0+0 _ s'tE') 
J detO J 



2[det0] 

x[det(0t0)] 2 0( 



detO 



-det0<5(l- , , 
2 v detO 



y/det[OW] 
det0+ 



).F[0]5(E'-£'t) 



) no\ 



J PS' 5(0+0 - E' t E')<5(E' - E't) 



-detO 5(0) T[0] 



J VY! 8{0 ] - E'+£')5(E' - E't) 



(A3) 



In the last step, we have used the property detO = det0t. 
Taking T[0] to be 

--^-=det0 y VY! 5(0+0 - E'+E'^E' - E't), (A4) 

eq. (A3) becomes 

VU 8{U^U - l)5(det[7- 1) 



= const, 



(A5) 



which is of the form of eq.(|18|) in the text. 

Next we look at the meaning of the variable E' in (A4) . 
The constraints on £' in (A4) are 



E't = £', E' 2 = 0+0. 



(A6) 



On the other hand, eqs.(|l3"|), ( |l4| ) and ( p0| ) in the text show 
that the cr-field is constrained as 



o-t = cr ; (ft+crfi) 2 = O^O. 



(A7) 



Comparing (A6) with (A7), we find £' ~ n+ufi. We know 
that the definition of a is not unique. It is up to a hidden 
symmetry transformation a — > h) ah. Therefore E' can be 
regarded as an equivalent definition of cr, and thus (A4) can 
be written as 



j!— = dct0 y Va 5(0+0 - rfo)6((T - o-t), 



T[0] 



which is just eq.(|T^) in the text. 



(A8) 
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